Entangling atoms in bad cavities 
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We propose a method to produce entangled spin squeezed states of a large number of atoms 
inside an optical cavity. By illuminating the atoms with bichromatic light, the coupling to the 
cavity induces pairwise exchange of excitations which entangles the atoms. Unlike most proposals 
for entangling atoms by cavity QED, our proposal does not require the strong coupling regime 
g 2 /kT 3> 1, where g is the atom cavity coupling strength, k is the cavity decay rate, and F is the 
decay rate of the atoms. In this work the important parameter is Ng 2 /hF, where N is the number 
of atoms, and our proposal permits the production of entanglement in bad cavities as long as they 
contain a large number of atoms. 

PACS numbers: 03.65.Ud,03.67.-a,42.50.-p 



I. INTRODUCTION 

To obtain a large coherent coupling of individual quan- 
tum systems while at the same time maintaining a low 
decoherence rate is the main challenge in the experimen- 
tal exploration of entanglement. An example of this is 
cavity QED which was one of the first proposals for the 
construction of a quantum computer and creation of en- 
tanglement of atoms [jf], but where the experimental 
progress has been hampered by decoherence caused by 
cavity decay and spontaneous emission from the atoms. 
To overcome these problems experiments have resorted 
to very small optical cavities where the small cavity vol- 
ume increases the interaction strength 0, 0] , or Rydberg 
atoms in superconducting microwave cavities, where the 
decoherence rates are low ||, |(| . 

Most cavity QED schemes which have been proposed 
so far require the cavity to be in a strong coupling regime 
g 2 /kT ^> 1 , where g is the atom cavity coupling strength, 
k is the cavity decay rate, and T is the decay rate of the 
atoms, and this limit is very hard to achieve experimen- 
tally. In this paper we propose a scheme where the cre- 
ated entanglement depends on the parameter Ng 2 /nT, 
where N is the number of atoms. With this scheme it is 
in principle possible to produce entanglement in any cav- 
ity, but in practise the entanglement becomes unmeasur- 
able if Ng 2 /nT < f . On the other hand, if Ng 2 /nT > 1 
a measurable entanglement is produced. Compared to 
the requirement g 2 / 'kT 3> f the present approach thus 
allows a substantially reduction in the requirements for 
the cavity if a large number of atoms is used. We note 
that the requirement Ng 2 /nT > 1 is equivalent to the 
criterion for optical bistability and squeezing in cavity 
QED as studied experimentally in Refs. |[ |). 

We propose to produce so-called spin squeezed states 
fiof . The collective properties of N two level atoms 
are conveniently described by pseudo angular momen- 
tum operators defined by J z = J2k(\ a )k( a \ ~ \b}k(b\)/2 
and J + = J2k \ a )k(b\, where the sum is over the indi- 
vidual atoms, and where \a) and \b) are the two internal 



states of the atoms. The state where all atoms are in the 
a state is an eigenstate of the J z operator with eigenvalue 
N/2. If the J x operator is measured in this state, the re- 
sult will fluctuate around the mean value of zero with 
a variance N/4. By entangling the atoms it is possible 
to maintain a large value of the mean spin (J 2 ) w N/2 
while considerably reducing the noise in a spin compo- 
nent Jg = cos(9)J x + sm(6)J y perpendicular to the mean 
spin. A state with this property is called a spin squeezed 
state. 

The experimental generation of spin squeezed states 
has potential applications in high precision spectroscopy 
and atomic clocks. For spectroscopy on a collection of 
two-level atoms Wineland et al. |jff| have shown that 
the possible gain in precision by using a spin squeezed 
state is given by the quantity 



e 2 =min(^| 



(1) 



It has also been shown that £ 2 < f indicates that the 
atom are in an entangled state jl^, [l3) , and here we shall 
use £ 2 to characterize the entanglement of the atoms. 

Several schemes for the production of spin squeezed 
states have been proposed, and recently the first weakly 
squeezed states have been produced experimentally by 
absorption of squeezed light HL by QND-detection [jf5| , 
and by collisional interaction |jl6|. The possibility to pro- 
duce spin squeezed states by having a large number of 
atoms in a bad cavity was already proposed in Ref . p7| , 
but the method proposed here is more efficient; with the 
same cavity parameters our proposal enables a stronger 
squeezing of the spin. During the preparation of this 
work we became aware that a scheme very similar to ours 
has recently been proposed |f8| . Compared to that work 
we apply a simpler level scheme in a different regime, 
but the fundamental ideas and the results are very sim- 
ilar. Also the results reached in Ref. |l!| are similar to 
ours, but the mechanism employed in that paper is very 
different from what we propose here. 

The paper is organized as follows. In Sec. || we present 
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our scheme in the ideal case where no dissipation takes 
place. In Sec. Ill we analyse the scheme in the more 
realistic situation where cavity decay and decay of the 
excited atomic states affect the preparation of the en- 
tangled state. In Sec. ^ we summarize our proposal, 
and we discuss how some states are more robust than 
others against dissipation and loss. The derivation of the 
evolution in the presence of dissipation uses the standard 
method of adiabatic elimination and produces quite com- 
plicated equations. These expressions are not essential 
for the understanding of the functioning of the proposal 
and we have put the technical derivation in the Appendix. 



II. IDEAL CASE 

The energy levels of the atoms and the laser couplings 
are depicted in Fig. [I] (a). We consider a A type three 
level atom with two stable ground states \a) and \b) with 
an energy difference uj ab and an excited state |e) with en- 
ergy difference uj a e to the ground states \a) (h = 1). The 
state | a) is coupled to the excited state |e) by a laser with 
a resonant Rabi frequency f2i and a frequency uj\ which is 
detuned from the excited state. Similarly the state \b) is 
also coupled to the excited state by another detuned laser 
with resonant Rabi frequency Q 2 and frequency w 2 . The 
two frequencies of the lasers are chosen such that their 
difference is exactly twice the energy splitting between 
the two ground states lu± — uji — 2uj ab . With this choice 
of frequencies all transitions involving only a single atom 
are off- resonant, but a transition which transfers pairs of 
atoms from state \a) to |6) is resonant. A similar choice of 
resonance conditions has also been proposed for trapped 
ions |2^, |TJ, and recently this scheme has allowed the 
first experimental production of four particle entangled 
states Q. To produce the pairwise excitations of the 
atoms, it is not sufficient that the process is resonant; 
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FIG. 1: Energy levels and couplings, (a) The energy levels 
of an atom and the couplings induced by the lasers and the 
cavity. The laser and cavity couplings permits the pairwise 
transfer of atoms from state \a) to |6) as indicated in part (b). 
The dashed lines indicate transitions induced by the cavity 
and the full lines are the laser couplings. 



it is also necessary that there exists a physical mecha- 
nism which enables an interaction between the atoms. 
In |2^, this was done by the Coulomb interaction be- 
tween the ions. Here we assume that the quantized field 
in an optical cavity couples both the states \a) and \b) to 
the excited state |e) with coupling constants g a and g b 
respectively. With this coupling to the cavity there exists 
a transition path for the pairwise transition, as shown in 
Fig. [I] (b) , and the matrix element for the transition be- 
comes non-zero. In the remainder of this section we show 
that this coupling leads to a spin squeezed state if we ap- 
ply the coupling to a state where all atoms are initially 
in the a state. 

If we assume all fields to be propagating in the same 
direction, the experimental situation is described by the 
Hamiltonian 



N 



H = u c!c+ y^cj ge |e) fe (e| + uj ab \b) k (b\ + H iuttk 
fe=i 

Y e-<Wlt +5-*)|e>*< | (2) 



le- iU2t + g b c)\e) k (b\+H.C., 



int.k 



where c and luq denote the annihilation operator and fre- 
quency of the relevant cavity mode. 

If we are in a regime where the laser power is suffi- 
ciently weak that we do not transfer any population to 
the excited atomic state we may adiabatically eliminate 
this state and obtain an effective Hamiltonian for the cou- 
pled state of the ground states and the cavity J2^] . As- 
suming further that the lasers are also sufficiently weak 
that we do not create a significant photon excitation in 
the cavity we may also adiabatically eliminate the cavity 
field and we are left with an effective Hamiltonian for the 
atoms il 



H 
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(3) 



where we have omitted some unimportant energy shifts, 
and we have introduced the detunings from the excited 
state Ai = ujae-^i and A 2 = w ae -w a i)-w2 = Ai+cj Q 6, 
and the detuning from the cavity mode 5 = u>i — uj ab — 
lu = uo 2 + ^ab — ^o- These detunings are also defined in 
Fig. ^. The angular momentum operators are defined as 
the angular momentum operators in Sec. I. The origin of 
each of the terms in this Hamiltonian can be understood 
from processes like the one shown in Fig. [I] (b) which 
gives the term with J+J+, i.e., a double Raman process 
which takes two atoms from \a) to \b) by absorption by 
f2i, emission into the cavity by <?&, reabsorption of the 
cavity photon by g a , and emission by 17 2 . 

If we choose the strength of the two Raman processes 
to be identical Q^/At = fi 2 .g*/A 2 = fig* /A, Eq. 
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(||) reduces to the simpler form -ffidcai = where 
X = |£l| 2 |.g| 2 /A 2 <5. The squeezing arising from this 
Hamiltonian can be calculated analytically jlTJ . Starting 
from an initial state where all atoms are in the a state 
and propagating with this Hamiltonian, squeezing by a 
factor of £ 2 ps N~ 2 / 3 is produced (in the limit N ~> 1), 
and this is a significant noise reduction if a large number 
of atoms is present. In the following section we show that 
a significant squeezing is produced even in the presence 
of dissipation. 

The time it takes to produce a spin squeezed state 
of many atoms is very short. To squeeze the spin by 
a constant factor, a constant number of atoms has to 
be transfered into the b state. With increasing N a de- 
creasing fraction of the atoms has to be transfered, and 
thus a shorter time (scaling as 1 /N) is necessary to make 
the squeezing. The different decoherence mechanisms 
therefore have less time to affect the preparation of the 
squeezed states, and as we show below, this reduces the 
experimental requirement for the production of squeezed 
states. 



III. ANALYSIS INCLUDING DISSIPATION 
AND NOISE 

The main purpose of this paper is to demonstrate that 
it is possible to use a cavity to entangle atoms even in 
situations where substantial dissipation is present. In 
this section we analyse the performance of our proposal 
in the presence of the two main decoherence mechanisms: 
spontaneous emission and cavity decay. 

Before making a quantitative analysis of the effect of 
dissipation we first make a few simple estimates of the 
decoherence. For simplicity we shall here assume that 
the lasers have approximately the same Rabi frequen- 
cies (fi; ~ f2) and detunings (A; ~ A) and also that 
the cavity couplings are similar gi ~ g. The number 
of spontaneously emitted photons is estimated to be ap- 
proximately iV r ~ NTt\n\ 2 /A 2 , where T is the total de- 
cay rate. The time required to produce squeezing by a 
constant factor is given by t ~ A 2 5/(Ng 2 H, 2 ), and by in- 
serting this expression we find that the total number of 
decayed atoms is 



(4) 



The number of photons decaying out of the cavity during 
the same time is estimated to be N K ~ Nntfl 2 g 2 / (A 2 S 2 ) 
which reduces to 



K 

a 



(5) 



when the expression for the time is inserted. 

Because the spin squeezed state are only weakly entan- 
gled they are quite insensitive to the spontaneous emis- 
sion of the atoms. If we assume that we are near the 
initial state where J z ~ N/2, J x and J y may be replaced 



by the canonical conjugate position x = J x y2 /N and 
momentum p = J y y/2/N operators of a harmonic oscil- 
lator. The relaxation rate for the harmonic oscillator is 
then the same as the relaxation rate for a single atom, 
and from the well known properties of squeezing of har- 
monic oscillators we find that the squeezing is not com- 
pletely degraded as long as the number of decayed atoms 
is much less than the total number of atoms. The de- 
cay of photons out of the cavity is more severe than the 
decay of a single atom. Because the cavity couples to a 
collective degree of freedom, the decoherence of the cav- 
ity will also affect the collective degree of freedom. If 
|fii<7b|/Ai = |fi2<7<,|/A2 we estimate that the first decay 
of a photon out of the cavity increases the variance in all 
directions perpendicular to the mean spin by a factor of 
three. To obtain a large squeezing we therefore require 
that at most a few photons are scattered out of the cav- 
ity. From the expression in Eqs. (Q) and (||) we see that 
we can fulfill both N-p -C N and N K < 1 if the cavity pa- 
rameters fulfill Ng 2 3> kT and we thus expect to be able 
to produce substantial squeezing in this regime. This is 
confirmed by our more accurate treatment of dissipation 
below. 

Dissipation is described by the master equation for the 
density matrix p 



dt 



p = -i 



H, p] + ^ {^.dmpdX, - d) m d m p - pd^dr^j , 



(6) 



where d m are relaxation operators. We assume that the 
separation of the atoms in the cavity is much larger than 
the wavelength of the spontaneously emitted photons. 
In this limit the decay of the atoms is uncorrelated and 
can be described by independent relaxation operators for 
each atom. The excited state |e) is assumed to have 
three independent decay channels: it may decay to the 
two lower states in the A-system |o) and \b) with decay 
rates j a and 71, respectively, and it may decay to some 
other state |o) with a decay rate j . The total effect of 
the spontaneous emission is described by 3iV relaxation 
operators 



d a , k = \^u\a)k(e\ 

db, k = Vlb~\b)k(e\ 
d Q ,k = y/%\o) k (e\, 



(7) 



where k — 1, N. To describe the decay of the cavity 
with a rate k we introduce a relaxation operator 



d r = 



(8) 



To derive the evolution of the spin squeezing we first 
adiabatically eliminate the excited state assuming 



< A 2 + — 



5, k' <C UJab 



(9) 
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FIG. 2: Evolution of squeezing for N = 10 atoms in a bad 
cavity (full line). The parameters used in the simulation are 
r = k = WOg, 7a = 7fc = To, Hi = «2 = 10 4 ' 



9a = = g, 



Ai = 10 J g, 5 = 5 • 10 g, and tj a 5 = 10 g corresponding to 
g 2 /kF = 10~ 4 but Ng 2 /kT = 10 2 . For comparison we also 
show the evolution with the same parameters but without 
dissipation, F = k = (dashed line). 



where I — 1,2, and T is the total decay rate of the ex- 
cited state r = 7 Q + 7b + j c . k' is an effective decay 
rate of the cavity, which is slightly larger than k due to 
the scattering of cavity photons by the atoms, cf. Eq. 
(A4). We then adiabatically eliminate the cavity from 
the equations assuming 



(10) 



Fig. U we show the evolution of squeezing in a situation 
where g 2 /Tn < 1 but Ng 2 /Tn > 1. With the realis- 
tic parameters used in the figure we are able to produce 
squeezing by approximately an order of magnitude after 
a very short interaction time. The time required to make 
the squeezing in the figure is less than 1/is if we chose a 
realistic cavity coupling parameter g = (27r)100kHz. 

In principle the proposed scheme can be used to pro- 
duce entanglement in any cavity. The field only cou- 
ples to the ground state |b) that is initially unpopulated, 
so that a large value of H2 increases the rate of the co- 
herent transfer of atoms from state a to & but does not 
affect the initial decoherence rate. The amount of squeez- 
ing, however, depends on the cavity parameters. To in- 
vestigate the obtainable squeezing we have performed a 
numerical optimization of the coupling strengths and de- 
tunings for a number of different cavity parameters. With 
fixed values of the dissipation rates and energy difference 
ujab, we vary 02/^1, S, and Ai and search for the values 
which give the minimal £ 2 (because all terms in M in- 
volve the square of the field strength, the minimum only 
depends on the ratio between the two fields). In the limit 
N ^> 1 and for fixed ratios between the decay rates 7; 
(/ = a, b, o) the results of the optimization indicate that 
the optimal squeezing parameter £^ in is only a function of 
the parameter Ng 2 /Tn. In Fig. ||we show £^ in for differ- 
ent cavity parameters. As expected from our simple es- 
timates, the figure confirms that strong squeezing can be 
produced in the limit Ng 2 /Tk 1. In the calculations 
we have assumed g a = gb = g and j a = jb = ~f = T/3, 
and with these values the opt imal squ eezing is approxi- 
mately given by £^ in = 0.7/^/Ng 2 /TK for Ng 2 /Tn > 1. 
This is indicated by the dashed line in the figure. The 



The equations resulting from the adiabatic elimination 
are quite complicated and we leave the derivation to Ap- 
pendix [A]. If we assume that the initial state is almost 
unaffected by the interaction so that J z w N/2, the ma- 
trix elements which are quadratic in the angular momen- 
tum operators do not couple to higher order terms and 
we obtain closed equations for the expectation values 



d_ 

dt 



(N a + N b ) 
(J+J+) 
(J-J-) 
(J+J-) 
(J-J+) 



M 



(N a + N b ) 
(J+J+) 
(J-J-) 
(J+J-) 
(J-J+) 



(11) 



where Ni is the number operator for atoms of type I = a,b 
{{N a + Nb) is n °t a conserved quantity because of the 
decay to the state \o)). The precise form of t he matrix 
M is given by the expressions in Eqs. (A5l A9). 

Due to the complicated structure of the matrix M it 
is difficult to describe the evolution of the squeezing an- 
alytically, but it is straightforward to find the evolution 
numerically, where the solution at a given time t may by 
found by taking the exponential of the matrix Mt. In 




0.1 



0.01 



NgVxT 

FIG. 3: Minimum squeezing parameter obtained by a nu- 
merical optimization. The points (+) are the results of the 
minimalization and the dashed line 0.7 / \/ Ng 2 / kT approx- 
imates for Ng 2 /kT > 1. In the calculation we have 
assumed N = 10 6 , g a = gb = g, la = 7f> = To, and 
bJab = 10 . The same minimum is obtained for three different 
ratios k/T = 10 -2 , 1, and 10 2 , and the results are independent 
of N in the limit N > 1. 
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same behaviour but with a slightly different constant has 
also been found in |51| . The obtained results only change 
slightly if we vary the ratio between the coupling con- 
stants or between the decay rates. 

The results of the numerical simulations agree very well 
with the behaviour expected from the simple estimates. 
Our assumptions of weak excitation of the atoms and 
of the cavity held mode imply the necessary condition 
A/ 3> r. It turns out, however, that the detuning S 
from the intermediate state in Fig.l (b) with one cavity 
photon excited, does not need to be large. If only the 
coupling is weak enough VA^iW Ai | <§; k, and if the 
process that absorbs cavity photons is stronger than the 
one producing them, \Q, 2 g a \/(A% + T 2 /4) > \n 1 g b \/(Af + 
r 2 /4), the photon excited state can be eliminated, and 
we hnd good squeezing for all values of S. The minimum 
value of £ 2 is found for S = 0. 

Finally, lets us briefly comment on a few experimental 
aspect of our proposal. The analysis above shows that 
our proposal is robust against the spontaneous emission 
caused by the laser coupling and the decay of photons 
out of the cavity. The coupling of the lasers, however, 
introduces other possible decoherence mechanisms if we 
are not able to control the lasers with a high enough 
accuracy. In the treatment so far we have ignored the 
AC-Stark shifts caused by the lasers because they can be 
compensated by a small change in the frequencies. But 
the magnitude of the AC-Stark depends on the power of 
the lasers so that fluctuations in the power has a detri- 
mental effect on the squeezing. To suppress this effect we 
propose to adjust the relative strengths of the two fields 
so that IftipA^A 2 + T 2 /4) = |ft 2 | 2 A 2 /(A 2 + T 2 /4). 
With this choice the AC-Stark shifts of the two ground 
states become identical and have no effect on the internal 
state preparation. The problem of stabilizing the power 
can thus be reduced to the problem of stabilizing the rel- 
ative frequency and intensity of the two fields which is 
much easier experimentally if the two fields are derived 
from the same source. 

The efforts to entangle atoms through cavity QED 
have so far concentrated on the strong coupling regime 
g 2 /nT ^> 1. To achieve this limit it has been desirable 
to use very small standing wave cavities where the cou- 
pling constant varies sinusoidally along the cavity axis 
with a period of half the optical wavelength. A con- 
trolled evolution in these cavities therefore requires that 
the atoms are localized in regions smaller than the opti- 
cal wavelength. Since our proposal puts much less strin- 
gent requirement on the cavity parameters it should not 
be necessary to use such small cavities and it could for 
instance be implemented with a ring cavity. Then, the 
magnitude of the coupling constant does not depend on 
the position of the atom along the cavity axis, and if the 
classical fields are co-propagating with the cavity field 
it is no longer necessary to localize the atoms within a 
wavelength. The atoms only need to be confined within 
the waist of the cavity mode, and this can be done with 
cold atoms trapped in a far detuned optical dipole trap 



or optical lattice, or even with atoms in a glass cell at 
room temperature. 



IV. CONCLUSION 

We have shown that it is possible to observe signifi- 
cant spin squeezing of atoms coupled to the field mode 
in a lossy cavity. The loss of quantum correlations be- 
tween the particles which is caused by atomic decay and 
cavity loss is balanced by the strong non-linear coupling 
achievable in the limit of very many atoms. It has been 
shown that spin squeezing implies entanglement, i.e., a 
separable state cannot lead to values of £ 2 smaller than 
unity. We have thus created entangled states which are 
fairly robust against dissipation and loss. 

There is no precise quantitative measure for the entan- 
glement of a large collection of particles, but a natural 
qualitative measure is to consider the possible gain, e.g., 
in spectroscopic resolution, that the entangled states of- 
fers with respect to a disentangled state. By binding 
the TV atoms together in N/P maximally entangled P- 
particle states (\aaa...a) + \bbb...b))/y/2, one obtains a 
spectroscopic resolution corresponding to a state with 
£ 2 = 1/P j24j]. Hence in terms of spectroscopic reso- 
lution the spin squeezed states are as powerful as if the 
atoms had been divided into groups of maximally entan- 
gled states of l/£ 2 particles. 

From a practical perspective, however, the squeezed 
states offer a significant advantage compared to a collec- 
tion of highly entangled states. In spin squeezed states 
the relevant observables are collective operators involv- 
ing all the atoms. There is no need to address the 
atoms individually, and the manipulation and detection 
of the squeezing can therefore be achieved by lasers ad- 
dressing all atoms collectively. Furthermore the spin 
squeezed states are also easier to produce: In an ideal 
spin squeezed state the one-particle density matrix p\ 
is very close to the initial pure state projection opera- 
tor pi = \a)(a\ + 0(1 /N), and the state of each atom 
is thus almost disentangled from the state of the other 
atoms. This means that we only need to perturb the ini- 
tial state slightly to turn it into a squeezed state. This 
is a significant advantage in any experimental attempt 
to produce entanglement because (a) the states are more 
robust against decoherence than more highly entangled 
states and (b) the interaction time required to make the 
desired state is much shorter than for the highly entan- 
gled states. The large number of atoms increases the 
decoherence rates, i.e., more photons are scattered, but 
our calculations show that the two advantages (a) and (b) 
outweigh the increased decoherence rate and enable the 
construction of entangled states in situations where the 
experimental capabilities do not permit the construction 
of entangled states of a few atoms. 

In a broader context our proposal fits into the field 
of ensemble quantum information processing, where the 
quantum information is encoded into the collective de- 
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grees of freedom of a collection of atoms. A number of 
papers [^5[ ^(| have proposed schemes for the pro- 
cessing of information encoded in such a way and the 
first experimental implementation of these concept has 
recently been reported [p8| . By combining these ideas 
with the present work one may for example imagine a 
quantum computer with several separate and individu- 
ally addressable atomic clouds which communicate in a 
controlled manner via cavity modes. 
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APPENDIX A: DERIVING THE EQUATIONS OF 
MOTION 



In this section we derive the equations describing the 
time evolution of squeezing in the presence of dissipation. 
We first consider only the Hamiltonian describing a single 
atom and we adiabatically eliminate the excited state of 
the atom by assuming that the population of that state 
is negligible. In this approximation the equations for the 
ground state density matrix elements are equivalent to 
the evolution by a Hamiltonian 



H = - 
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and six relaxation operators dk,i {k = a, 6, o and / = 1,2) 
describing the combined excitation with a detuning A; 



and decay to a state \k), e.g., 



d a ,i 



Ha 



Hi 
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(A2) 



To derive these result we have assumed that 5 <C A 1; A2, 
and we have used the second relation in Eq. (|J) to neglect 
processes which creates photons without a change in the 
atomic state. 

The first two lines in Eq. (Al) represent AC-Stark 
shifts of the ground states. The first part of the shifts 
containing the classical fields fii and f2 2 can be compen- 
sated if we make a change in the frequency of the fields. 
The second part containing the quantum field c is much 
smaller than the first and by inserting the approximate 
time and Eq. (A3) below, we find that this term gives a 



negligible phase shift if g 2 /5A <C 1 and we shall neglect 
these terms. 

We then adiabatically eliminate the cavity in the 
Heisenberg picture. Setting d(ce iSt )/dt = we obtain 



ASt 




5 + if 



noise. 



(A3) 



where the noise ensures the commutation relation of the 
operator. Here we have introduced an effective decay rate 
for the cavity 
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(A4) 



which takes into account that the cav ity photons may be 
scattered by the atoms. In Eq. (A4) we have assumed 
that essentially all atoms remain in the a state. 

The adiabatic elimination of the cavity requires that 
the atoms are disentangled from the cavity, i.e., that 



c) <C 1. At t = this reduces to Eq. (10), and at 



later times our numerical simulation indicate that (c*c) 
typically changes slowly, so that the condition is fulfilled 
if it is fulfilled at t = 0. 



Finally, we insert Eq. (A3) into the time derivatives of 



the angular momentum operators, and by assuming that 
the initial state only changes slightly so that J z sa N/2, 
we obtain the following expressions 
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